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We study the effects of hybridization between the two electron pockets in Fe-based superconduc¬ 
tors with s-wave gap with accidental nodes. We argue that hybridization reconstructs the Fermi 
surfaces and also induces an additional inter-pocket pairing component. We analyze how these two 
effects modify the gap structure by tracing the position of the nodal points of the energy dispersions 
in the superconducting state. We find three possible outcomes. In the first, the nodes simply shift 
their positions in the Brilluoin zone; in the second, the nodes merge and disappear, in which case 
the gap function has either equal or opposite signs on the electron pockets; in the third, a new set 
of nodal points emerges, doubling the original number of nodes. 


I. INTRODUCTION 


The iron pnictides and chalcogens have been the sub¬ 
ject of intense study since 2008, when it was discovered 
that they are superconducting at relatively high critical 
temperature^. Understanding their gap structure is an 
important step toward identifying the mechanism respon¬ 
sible for superconductivity in these materials. Although 
the multi-band nature of Fe pnictides/chalcogenides al¬ 
lows for many different gap structures, the prevailing sce¬ 
nario is that the pairing occurs between electrons on the 
same Fermi surface (FS) and the superconducting gap 
function has s~^ symmetry, i.e., the gap changes sign 
between hole and electron pockets. There is experimen¬ 
tal evidence that in some members of the family, like 
BaFe 2 (Asi_a,Pa,) 2 ,^ LaOFeP,!^ and LiFeP,® the gap has 
nodes, likely on the electron pockets. 

Previous studies of the gap structure were mostly re¬ 
stricted to an Fe-only approach, in which a generic model 
of the band structure consists of two nearly circular hole 
pockets centered at (0, 0) and two elliptical electron pock¬ 
ets centered at (tt, 0) and (0, tt) in the first Brillouin zone 
(BZ) (see Figure]^. In some systems there exists, at least 
for some kz, a third hole pocket, centered at (7r,7r). The 
s~^ superconductivity is believed to be chiefly caused 
by a magnetically enhanced interaction between hole and 
electron pocket^. The nodes on the two electron pock¬ 
ets come about because by symmetry the s-wave gap 
on these pockets has the form A(I ± a cos 20k) (plus 
higher harmonicsP and if a > I, the gap vanishes at 
cos2Ok = ±l/a0^. 

However, this Fe-only scenario is incomplete because 
the electron hopping between Fe atoms predominantly 
occurs via pnictogen or chalcogen sites, half of which are 
located above and half below each Fe layer in a checker¬ 
board patteriP. As a result, the actual symmetry is lower 
than that of the Fe-only lattice, and the correct unit cell 
contains two Fe atoms. The non-equivalence of hopping 
from above and from below an iron layer causes the elec¬ 
tron pockets to hybridize. In this paper we will study the 
effect of this hybridization on the gap structure. 

In order to incorporate this effect into models with 
one Fe atoms, one has to include additional terms in the 



Figure I: Unfolded Brillouin zone with one Fe atom per unit 
cell (with no hybridization). Holes and electron pockets are 
labeled by h and e, respectively. The crosses on the electron 
pockets indicate where the nodal lines of the gap function 
intersect the Fermi surface. 


Hamiltonian with excess momentum Q = (7r,7r). This 
does not actually violate conservation of momentum be¬ 
cause this vector folds into a reciprocal lattice vector in 
the actual BZ with two Fe atoms per unit cell. The mo¬ 
mentum Q connects the two electron FSs and also the 
hole pocket centered at (7r,7r) with the other two hole 
pockets. Our primary goal will be to study how the ac¬ 
cidental nodes on the electron pockets evolve once we 
include hybridization. Therefore, we focus on the effect 
of hybridization on the electron pockets. 

The hybridization gives rise to two effects. First, hop¬ 
ping via a pnictogen/chalcogen sites gives rise to an addi¬ 
tional quadratic term in the Hamiltonian for two electron 
pockets 


^a = E 

k 





( 1 ) 


where c and d are operators for electrons near each of 
the two electron FSs (we discuss the form of Xk in the 
next section) and the sum over repeated spin indices is 
implied. This cross-term mixes the two electron pock¬ 
ets and reconstructs the electron FSs. Second, there ap- 
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pear new four-fermion interaction terms in which incom¬ 
ing and outgoing momenta differ by Q. In the super¬ 
conducting state, in which we are interested, two out 
of four fermions can be put into the condensate and 
the four-fermion terms with excess momentum Q reduce 
to quadratic terms with prefactors proportional to the 
superconducting gap. These new terms describe inter¬ 
pocket pairing between fermions from two different elec¬ 
tron pockets: 
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^kQ!^-k-Q/3 


+ di 


ka^-k-Q/3 


la: 


cxf3 


+ H.C. 

( 2 ) 


In other words, due to hybridization, the non¬ 
zero intra-pocket pairing condensates 

(^k+Qa^-k-Q/ 3 ) iiiduce inter-pocket pairing between the 
two electron pockets. 

In this paper we study how the additional terms 1-Lx 
and affect the gap structure when nodes are present 
on the electron pockets. The effect of the hopping Ak 
term alone has been studied befor^^^, but not its inter¬ 
play with the pairing term. We find that the hopping 
and inter-pocket pairing terms generally pull the nodal 
points in opposite direction. If the Ak term dominates 
and reaches a certain threshold value, the nodes merge 
and disappear at particular symmetry points, and the 
gap acquires a uniform and equal phase on the two elec¬ 
tron pockets (opposite to the phase on the hole pockets). 
In contrast, when f3k dominates and reaches a thresh¬ 
old value, the nodes merge and disappear at a different 
set of symmetry points and the phase of the supercon¬ 
ducting order parameter becomes opposite on the hy¬ 
bridized electron pockets. This is the same gap structure 
that was recently found in the analysis of pairing in the 
orbital formalisnJ^ and dubbed orbital anti-phase s~^ 
state. The state with opposite signs of the gaps on the 
two electron pockets has also been found in the analy¬ 
sis of possible superconducting states in LiFeAs, albeit 
for a different reasorP^. As an interesting peculiarity, 
we found that for elliptical pockets nodes disappear in a 
rather non-trivial way - first new nodes appear and the 
number of nodes doubles, and then the new and already 
existing nodes merge and disappear. Such behavior has 
not been found before in the studies of multi-band super¬ 
conductors, as far as we know. 

The hybridization between electron pockets, either due 
to Hx or to H/ 3 , has to be distinguished from the effect of 
the folding of the IFe BZ into the 2Fe BZ. Upon folding, 
fermionic momenta transform as = kx ^ ky^ ky = ky — 
kx^ and the two electron FSs, originally centered at (0, tt) 
and (7r,0), merge around (7r,7r) (see Fig. [^. The inner 
and outer FSs touch each other along kx = tt or ky = tt. 
The merging can be viewed as “reconstruction” of the 
two electron FSs into an inner FS with no nodes and 
an outer FS with 8 nodes, or vice versa. However, this 
“reconstruction” is just a rotation in momentum space 
and a re-branding. The kp remains the same and the 
location of the nodes remains at exactly the same angles 


Unfolded BZ Folded BZ 




Figure 2: Fermi surface folding, a) Unfolded Brillouin zone 
corresponding to a IFe unit cell, b) Folded Brillouin zone 
corresponding to a 2Fe unit cell. The electron Fermi surfaces 
overlap and are reconstructed into outer and inner parts. De¬ 
pending on their original location, the nodes of the gap func¬ 
tion are either all located on the outer surface, as shown in 
c), or on the inner surface, as shown in d). 


as without folding, only the reference axis rotates by 45*^. 
The hybridization is a different phenomenon - it actually 
reconstructs the original FSs into an inner and an outer 
FSs at new kp and creates a pairing component between 
them. As a consequence, the position of the nodal points 
shifts. 

The paper is organized as follows: In Section |TI| we dis¬ 
cuss our model. In Sectionjinj we consider, as a warm-up, 
the limiting case of circular electron pockets and analyze 
the two effects of hybridization first separately and then 
together and study their interplay. Section extends 
the analysis to the more general case of elliptical pock¬ 
ets. We summarize our findings in Section [V| 


II. THE MODEL 


We consider a two-dimensional multi-band model with 
hole pockets centered at (0,0) and elliptical electron 
pockets centered at (0, tt) and (tt, 0) in the unfolded BZ. 
In the normal state, the Hamiltonian describing the two 
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non-hybridized electron pockets is simply 

^0 = + ^k<^ka<^kQ! 

k 


( 3 ) 


We assume that the dominant interaction which leads 
to s~^ superconductivity is the repulsion between elec¬ 
tron and hole pockets, enhanced by (tt, tt) spin fluctua¬ 
tions (in the 2Fe BZ). The s-wave gap on the hole pockets 
is invariant under rotations by 7r/2, so it can be expanded 
in cos4n0 harmonics, where (j) is the angle along the 
hole pocket and n is an integer. On the electron pockets 
the expansion of the 5—wave gap contains cos 2n0 terms, 
where 0 is the angle along the electron pockets and the 
components with odd n = 1, 3,... have opposite signs on 
the two electron pockets, if we measure 0 from the same 
direction for both pocket^. These odd multiples of 20 
are allowed because the electron pockets transform into 
each other under a rotation by 7r/2. Numerical analysi^^ 
shows that the gap on the electron pockets can be well- 
approximated by the two first harmonics n = 0 and n = 
1, whose magnitudes are generally comparable to each 
other. Accordingly, we set Aei(^k) = A(1 — acos20\^), 
^e 2 (^k+Q) = A(1 + (acos26>k+Q). The corresponding 
term in the Hamiltonian is thus 


"Ha = (1 -yk) 


J J 

‘'ka^-k/3 


( 4 ) 


+ (1 + yk+Q)4+Qa^-k-Q/3 


^(7 


aj3 


+ H.C., 


where y\^ = a cos 26>k and ^k+Q = cos 26>k+Q. 

Our goal is to analyze how accidental nodes on electron 
pockets evolve with increasing hybridization. For this we 
assume from the beginning that \a\ > 1 in which case 
the gaps on electron pockets have nodes when ^k = 1 
(for c—fermions) and ^k+Q = for d—fermions. 

To simplify the presentation, we fold the IFe BZ into 
the 2Fe BZ and replace the momentum k + Q of d- 
fermions by k. The momenta k below are defined as 
a deviation from Q, which is the location of the electron 
pockets in the folded BZ. In order to preserve the cos 2dk 
form of the gap function, we define dk relative to the 
minor axis of the c pocket. 

In the normal state, the inclusion of hopping via pnic- 
togen/chalcogen atoms generates mixing between c and 
d fermions: 




( 5 ) 


A microscopic derivation of Ak show^^^^I^®^ that in 
1111 systems (in which the configuration of pnictide 
atoms around every Fe layer is the same), Ak vanishes 
along the lines = -^ky and has some weak kz depen¬ 
dence. In 122 structures (in which the “above/below” 
configuration of pnictogen/chalcogen atoms is inverted 
from one Fe layer to the other), Ak has minima but does 
not vanish along any direction. In the presence of a spin- 
orbit interaction Ak does not have zeros even in 1111 sys- 
tem^IHIIll. Because our goal is to understand the generic 


effect of the hybridizatio n betw een the two bands, we will 
treat Ak as a constant Earlier analysis of the ef¬ 

fect of Ak including its angular dependence (but without 
interplay with inter-pocket pairing) has shown that the 
results are qualitatively the same as for constant 

We next consider how hybridization affects the pairing 
terms. They can be subdivided into two types. Terms 
of the first type describe an interaction with excess mo¬ 
mentum Q between electron pockets and contain three 
fermionic operators from one pocket and one from the 
other pockefP^, e.g.. 



Terms of the second type contain an interaction with 
excess momentum Q involving two fermions from a hole 
pocket and two from different electron pockets, e.g.. 


+ hlfd 

p+q,/^) 


k,p,q 


, ( 7 ) 


where the operator hpa describes fermions near one of 
the hole pockets. The two types of terms are differ¬ 
ent, yet their effect on the s~^ superconducting state 
is the same - both induce an additional pairing inter¬ 
action between fermions belonging to different electron 
pockets. Indeed, in the superconducting state 
(^ka^-k/ 3 )’ and are all non-zero. Decoupling 

four-fermion terms in Q and 0 using these averages, 
we obtain anomalous quadratic terms involving c and d 
fermions: 


'^13 = 

k 


^ka^-k^ 


+ dL,c 


ka’^-k/? 


*^a/3 + H-C.. (8) 


The coupling Pk is proportional to the magnitude of the 
s~^ gap A and has some non-singular angular depen¬ 
dence which we can safely neglect. 

We assume without loss of generality that the parame¬ 
ters A and A are real and positive. The parameter /3 (x A 
is then also real, but its sign can be either positive or neg¬ 
ative. 

Below we consider various ratios of A/^d and two FS 
geometries. In each case we compute the quasi-particle 
dispersion in the superconducting state and determine 
the position of the nodal points. In all cases we find two 
different dispersions: One is gapped over the entire BZ, 
while the other contains nodal points in a subset of the 
parameter space. 
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III. CIRCULAR POCKETS 


B. Inter-pocket hopping only (A 7 ^ 0, /3 = 0) 


As a warm-up, consider the limiting case when the 
two electron pockets are identical and have full rotational 


symmetry, i.e., 


A. Inter-pocket pairing only (/3 7 ^ 0, A = 0) 


In this special case, we find, after straightforward di- 
agonalization of the quadratic form, that the two disper¬ 
sions in the superconducting state are 

(i^k)'=eU(A±^A2y2+^2^ , (9) 

where the expression in parenthesis repre¬ 
sents an effective gap function. Ai ^ = 0 , 

Et = ±v/g^ + A2(l + |yk|)^ and = 

ib-^/Ck + A2(l — |yk|)^- This corresponds to the 
gap structure in the folded BZ: one band has no nodes 
and the other band has 8 nodes (see Fig. |^. 

At /3 7 ^ 0, both dispersions evolve. The nodal points 
are still located on the FS (the locus of zero energy points 
in the normal state, given by = 0). The band with 
energy is shifted up at a non-zero [3 and its effective 
gap function is definitely nodeless. In contrast, the band 
with energy is shifted down and the positions of the 
8 nodes shift to 


a/A 2 - 

cos( 2 ^k) = ±^^-( 10 ) 

As P increases, the nodal points move toward the diag¬ 
onals kx = Eky, as shown in Fig. At a critical value 
/3c = A they meet along the BZ diagonals. At larger 
[3 > f3c the nodes disappear. At the same, because the 
nodes merge along the BZ diagonals, the sign of the gap 
on one FS is opposite to that on another FS. Such a 
gap structure has been obtained before in the analysis 
of possible gap configurations in LiFe As both in orbital 
formalism and in band formalisnffm^. In the orbital for¬ 
malism, such a state was termed “orbital antiphase” s~^ 
stat^m. 

An intuitive way to understand this behavior is the 
following: the gap function at the upper band is positive, 
while the one at the lower band has a roughly sinusoidal 
shape that crosses zero eight times, and its maxima occur 
at the diagonal directions kx = ky and kx = —ky. As \P\ 
increases, the gap function shifts downward and thus the 
nodes move towards the BZ diagonals, until \P\ reaches 
the critical value /3c. At this point, pairs of nodal points 
meet and annihilate. At larger |/3| > A the maxima of 
the gap function are located below zero, i.e., the gap is 
negative for all angles. This behavior is illustrated in Fig. 

SI 


This case has been studied before,!!^ and we briefly 
summarize the results for completeness. The hybridiza¬ 
tion of the electron pockets reconstructs the FSs. In or¬ 
der to study this effect we diagonalize the Hamiltonian 
3-Lq + Ex by introducing new quasi-particles a and b via 


Cka — 




<3kcK 


^kcn “1“ ^kcn 


( 11 ) 


After diagonalization, the non-pairing part of the Hamil¬ 
tonian takes the form 


k 


( 12 ) 


where the dispersions are = ^k + A and ^k ~ ~ 

A. The new FSs are concentric circles. The a pocket is 
smaller and the b pocket is larger than the unhybridized 
pockets. 

In order to study the superconducting state, it is con¬ 
venient to first rewrite Ha in terms of the new operators 
a and b as 




^ka^-k/3 + ^ka^-k/3 


la: 


a/3 


(13) 


■^Ayk 


®ka^-k/3 "ka“-k/3 


bt„a' 


+ H.c. 


Note that the inter-pocket pairing component A^k 
emerges. Diagonalizing the new Hamiltonian Hq + H'^, 
we find the two dispersions for the quasi-particles in the 
superconducting state given by 

=Ak±^, (14) 

where 

^k =Ck + A^(l -1- yl-) -f A^, (15) 

Bk =4 [(aA)2 + A2y2(A2 + a2)] , (16) 

The dispersion E^^ as defined in ( [l^ , is positive for all 
k even when evaluated at A = 0 , so hnas no locus of the 
nodal points. Both FS lines ^k = 3zX in the normal state 
are part of the other dispersion E^. Once A becomes 
non-zero, one can easily check that E^ is non-zero along 
the normal state FSs. However, on the original, non- 
reconstructed FS, 


E^ = v/A2 + A2 - A|yk| (17) 


This function contains 8 nodal points located at 

a/A^ -j- 

cos(20k) = 

As one increases A, the nodes stay on the unhybridized 
FSs, but move toward the kx or ky axes (whichever is 
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a) Original pockets 


b) Folded pockets c) Increasing inter-pocket pairing 



Figure 3: Gap structure for circular electron pockets with inter-pocket pairing. The crosses represent the location of the nodal 
points of the quasi-particle dispersions in the superconducting state. In a) and b) we show the unfolded and folded zones, 
respectively. In the folded zone, the pockets should overlap but we separate them for clarity. As we increase the inter-pocket 
pairing, keeping the conventional intra-pocket pairing fixed, the nodal points shift toward the diagonal lines kx = ^ky, as seen 
in c). If the inter-pocket pairing reaches a critical value, the nodes merge and disappear, resulting in electron Fermi surfaces 
with opposite signs of the gap function. 




CD 



0 -71/4 Till 371/4 71 

0 



Figure 4: Nodal gap function with inter-pocket pairing /3 eval¬ 
uated over the FS. As [3 increases, the function shifts down¬ 
ward and its zeroes move toward the angles n'Kj2^ where n 
is an integer, which correspond to the directions given by 
kx = 3zky. At /3 = A, pairs of zeroes merge at those angles. 
For yd > A, the function is negative and has no zeroes. 


closer), until A reaches a critical value Ac = — 1. 

At this value of A, pairs of nodal points merge and then 
disappear at larger A. We show this schematically in Fig. 

El 


The analysis of the signs of the gap requires some care. 
For A ^ A, the inter-pocket pairing term becomes irrele¬ 
vant as the two reconstructed FSs are far apart from each 
other. In this limit, the gap on each reconstructed pocket 
is given by the fist line in (13) and is just A for both pock¬ 
ets. In this limit, the sign of the gap is indeed the same 
on both FSs. At intermediate A, however, one cannot de¬ 
fine the phase of the gap function on the two FSs because 
in terms of the hybridized fermions the gap has contri¬ 


butions from both inter-pocket and inter-pocket conden- 
sates {ak,aa-k, 0 )i(^lp and {ak,ab-k, 0 )io-lp, respectively. 
Because the limiting behavior at large A is known, it is 
“natural” to define both finite gaps with the same sign, 
as in Eq. (13), but we caution that this is rigorously 
justified only in the limit of very large A. 


C. Inter-pocket pairing and hopping (yd 7 ^ 0, A 7 ^ 0) 


Once more, the first step is to diagonalize the Hamilto¬ 
nian Ho by introducing new pockets a and b exactly 
as in the case with yd = 0. We then rewrite the remaining 
Hamiltonian in terms of the a and b operators and obtain 


^A+/3 — 9 ^ /^)^ka^-k/3 






‘^ka^-k/3 + ^ka®tk/3 


ial^ + H.C. 


(19) 


Observe that the coefficient (3 appears only in the intra¬ 
pocket terms. We diagonalize the Hamiltonian + 
^A+/3 again find dispersions of the form = 

\/Ak ± where 

=ek + A2(l + 2/^) + A2+/32, (20) 

Bk =4 [(^A + A/3)2 + A2y/2(A2 + a2)] . (21) 

The dispersion is fully gapped, but has nodal 
points at momenta given by 


a = A/3/A, 


cos 20k = i 


v/A2 + A2-/32-A2/32/A2 
aA 


( 22 ) 

(23) 
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a) Unfolded b) Folded 



c) Increasing inter-pocket hopping 



Figure 5: Gap structure for circular pockets with inter-pocket hopping. In a) and b) we show the unfolded and folded gap 
structure in the absence of hybridization, respectively. Inter-pocket hopping reconstructs the FSs as shown in c). The dashed 
line represents the original pockets, while the solid lines are the new set of two pockets. As the hopping parameter increases, 
the inner FS shrinks and the outer one becomes larger. In the superconducting state, the nodal points (represented as crosses) 
lie on the unhybridized FSs and shift toward the kx and ky axis as the hopping increases. At a critical value of this parameter 
the nodes meet and merge in pairs. For greater values of the parameter they vanish and the sign of the gap function becomes 
the same on both FSs. The separation between the FSs is purely schematic and has been exaggerated for clarity. 


c) Increasing A (0<|P|<A) 


a) Unfolded b) Folded 



A: Inter-pocket hopping 
P: Inter-pocket pairing 



d) Increasing |p| (A>0) 



Large A 
|P|<A 


|P|>A 
Large A 


Figure 6: Gap structure for circular pockets with inter-pocket hopping and pairing interaction. The solid outer and inner circles 
represent the reconstructed FSs after hybridization. In the superconducting state, the nodal points (represented as crosses) lie 
on a circle (dashed line) but no longer on the unhybridized Fermi surface. As the hopping parameter A increases, the nodes 
move toward the kx and ky axis, where they can merge and disappear. In this scenario the gap function has the same sign 
on both electron pockets. Increasing the pairing strength [3 shifts the nodes toward the diagonal lines kx = 3zky instead. The 
nodes can merge and vanish, in which case the gap function acquires opposite signs on the outer and inner pockets. 


Note that the nodal points are now shifted from the un¬ 
hybridized FS. The direction of the shift depends on the 
sign of /3. If /3 > 0 (/3 < 0) the nodes appear between the 
original FSs (the ones before hybridization) and the outer 
(inner) hybridized FS. The nodes exist in the parameter 
range given by 

0 < +A2- A^/3VA^ < (24) 


The lower bound is reached when we keep A fixed and in¬ 
crease \/3\ towards critical Pc = In this case the nodes 
merge at the diagonals kx = ky or kx = —ky. The nodes 
disappear when \P\ > Pc and the intra-pocket gap com¬ 
ponents in the first line of Eq. (19) have different signs. 
In the limit of A ^ A, the inter-pocket gap component 
becomes irrelevant since the reconstructed FSs are far 
apart and in this sense the gap function has opposite 
signs on the two electron pockets. 

The upper boundary in (24) is reached when we set 
\P\ < A and increase A. In this situation the nodes merge 


along the kx and ky directions at a critical value of A given 
by 


Ar — A4 


' — 1)A^ + P‘^ 

A2 -^2 • 


(25) 


At A > Ac, the nodes disappear and the gap function 
has the same sign on each electron pocket, as can be 
clearly seen in the limit of A ^ A. Note that Ac grows 
with P^ i.e., the inter-pocket pairing allows the nodes to 
exist in a greater range of values of A. In this sense the 
pairing partially protects the nodes from disappearing 
due to hopping, as long as \P\ < A. The behavior of the 
nodes when both A and P are present is summarized in 
Fig. |6] 
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IV. ELLIPTICAL POCKETS 


Now we consider the more realistic case where the elec¬ 
tron pockets are elliptical. We will take the dispersions 
in the form 


?k = 
?k = 


-/i + 


A. 

2mi 


-/i + 


jA 

2m2 


2 m 2 ’ 


A_ 

2mi' 


(26) 

(27) 


It is convenient to rewrite the dispersions as i 

(5k cos 26>k, where the + sign corresponds to ^k = (Ck + 
^^)/2, and ^k = Without loss of 

generality, we will take (5k to be positive. 


A. Inter-pocket pairing only (/3 / 0, A = 0) 


The dispersions after diagonalizing the Hamiltonian 
are given by (^k i \/5k, where 

1 


^*^=2 1 


(C£)^ + (^£)^ + 2A^(l + y^) 

((^^)^-(e£)^+ 4 AV)' 

+4|/3iM(e£-e^)'+4A2) 


213^ 


(28) 


(29) 


Once again, is fully gapped but may contain 
nodes. Unlike the circular case, the nodes are not lo¬ 
cated on the original FSs, but at momenta |k| which are 
solutions of 

(^k + <^Ck) (^k^k - (aA^) -h (a/3^^k = 0. (30) 

Let the solutions to this equation be ^k = C ^k = 
The angular positions of the nodal points are given by 


cos^ 26>k = 


ms) 


(J 2 +c, 2 A 2 ) 2 ’ 


(31) 


where 


F(^, S) = ((5^ - a‘^A‘^) (C^ -h - A^) 

- AaA^^6. (32) 


Note that for each solution to Eq. (30) there exist 8 


nodal points in the dispersion. One may solve for ^ and 
S exactly but the solution is not very illuminating. It is 
more useful to solve for ^ in terms of 6 and analyze how 
the nodal points evolve when we change 6. Expressing f 
in terms of S we obtain 




6^ - a^A^ ± 


(a 2 A 2 + 


2a5 


(33) 


Substitution of these solutions into Eq. 


yields 


cos^ 26>k = 


52 - q;2^2 ^ ^(q;2A2 + j2)2 _ 4Q,2^2j2 


2a252 


( 34 ) 


Analyzing (34) we find new interesting physics. 


Namely, depending on the parameters, there may be 0, 
8, or 16 nodal points in the dispersion. When S < aA, 
there are either 8 or zero nodes, as one can immediately 
verify. At small (3^ there are 8 nodes. As |/3| increases, 
the nodes move toward the diagonals = Eky^ like in 
the circular case. At |/3| = A, pairs of nodes merge, 
and for \P\ > A they disappear. The outcome of the 
disappearance of the nodes is the effective supercon¬ 
ducting state with different signs of the gap on the inner 
and outer electron pockets, see Eig. 

When S > (aA, the evolution of the nodes is more 
interesting. At small /3^ there are again 8 nodes. As |yd| 
increases, the nodes shift towards diagonals but they do 
not reach k^ = Eky at |/3| = A. Instead, at this a 
new quadratic node appears in in each quadrant at the 
point where zone diagonals intersect the original ES. At 
I yd I > A, each quadratic node splits into two, one to the 
right and one to the left of a diagonal, and each new 
node moves toward the already existing nodes (see Eig. 
[^. Thus, there is a total of 16 nodal points. As |yd| 
continues increasing, the old and new nodes merge at a 
critical value \P\ = Pc given by 


Pc = (a^A^ + (5^) /{2\a\5). 


(35) 


The nodes disappear when |yd| exceeds this critical value 
and the end result of the evolution of the nodes is the 
same minus-plus gap on the inner and outer electron 
pockets. 

We verified this behavior numerically (Eig. [factually 
shows the results of numerical calculations). In all nu¬ 
merical examples here and below we have set fi = 10A 
and a = —1.5. Note in passing that while the relation¬ 
ships presented in this analysis are exact, one must keep 
in mind that in general 5 by itself depends on p. 


B. Inter-pocket hopping only (A / 0, ^ = 0) 

This case has been studied befor^^ and we present it 
here for completeness, using a somewhat different com¬ 
putation scheme. The first step is to diagonalize the 
Hamiltonian Ho + Ha by introducing new operators a 
and b such that 


Cka — ^k^ko! 4“ '^k^ko!? 
d\^a — '^k<^kcK 4“ '^k^kcK5 


where u\^ = cost/; and v\^ = —sin7/5, and this angle 
p) satisfies cos 27/; = (^g — Ck)/'\/(^k “ ^k)^ 

sin 27 /; = 2A/y'^(^^^-^^)M-4W. The new Hamiltonian 
is 


7^0 — [^k^^ka^ka 4 - ^k^ka^ka 

k 


(37) 


where A = + Ck) ± \/V + (^k -^k)V4- These 

new dispersions define the reconstructed ESs shown in 





















Increasing inter-pocket interaction 



New nodes emerge at 
Fermi surface crossing 


Figure 7: Gap structure for elliptical pockets with inter-pocket pairing. In a) we show the folded electron pockets overlapping. 
The crosses indicate the position of the nodal points, while the solid and dashed lines indicate opposite signs of the gap function. 
The evolution of the nodal points as obtained by a numerical calculation is shown in b) through f). First, b) shows the nodal 
points with no hybridization. If the ellipticity is below a threshold, increasing the inter-pocket pairing simply shifts the nodes 
toward the diagonal lines kx = :^ky as seen in c) and d). The nodes merge and disappear after reaching the diagonal lines and 
the gap structure becomes nodeless as shown in g). If this occurs, the gap function has opposite signs on the inner and outer 
edges of the FSs. Alternatively, if the ellipticity is higher than the threshold, the nodes do not reach the diagonal lines, but 
instead a node emerges at that symmetry point. As we further increase the inter-pocket pairing the node splits into two nodes 
which move toward the original nodes, as shown in e) and f). Eventually, the new and old nodes meet and merge, resulting 
again in the structure shown in g). 


Fig.H As the hopping parameter A increases, the outer 
FS (associated with b fermions) becomes larger, while the 
inner FS (associated with a fermions) shrinks. 

Now we consider the superconducting state. Rewriting 
the pairing part of the Hamiltonian in terms of the new 
operators we find 


In this particular case 



A. A 

'{^^? + ie^f+2A\l + yl) 

5 

(39) 

b,A_ 

4gb^-(e£)^+4AVkyk)' 




+ - eb" + 4a2) 

(40) 




I “ /k2/k)aL«tk/3 


+ ^(1 + /k2/k)&k„^Lk/3 


ta 


al3 


2 ^ 
k 


A^k^k 


<^ka^-k/3 


dL,c 


kcn^—k/3 


'a/3 


(38) 

H.C., 


where /k = cos 2-0 and g\^ = —sm2'ip. The diagonaliza- 
tion of this Hamiltonian again yields dispersions in the 
superconducting state in the form = Ak ± \/Sk. 


As usual, the dispersion is nodeless, but E^ has 
nodes at momenta which are the solutions of the equation 

{S\^ + Q^^k) (^k^^k - 0'A‘^) - aX^6\^ = 0. (41) 

Each solution to this equation defines a pair (^, J) and 
determines the radial position of the nodal point. The 
angular position is then given by 


where 


cos^ 26>k 


F{U) 

(52 + a2A2)2’ 


(42) 


F(e, 5) = {F - a^A^) - A" - A^) 

- 


( 43 ) 
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Increasing inter-pocket hopping 



Figure 8: Gap structure for elliptical pockets with inter-pocket hopping. In a) we show the gap structure in the absence of 
hybridization. Inter-pocket hopping reconstructs the FSs as shown in b). As the hopping parameter increases, the inner FS 
shrinks and the outer one becomes larger. In the superconducting state, the nodal points (represented as crosses) lie near, but 
not exactly on FSs and shift toward the kx and ky axis as the hopping increases, as shown in b) and c). At a critical value of 
the hopping the nodes meet and merge in pairs as in d). Finally, for greater values of the parameter they vanish and the sign 
of the gap function becomes the same on both FSs. Subfigures b) through e) are the result of numerical calculations. 


Like before, we solve for ^ in terms of 6. The solution 
is 


_ -52 + 

^-• (44) 


2a5 


When we substitute this solution into Eq. ( [4^ we find 
that the angular position of the nodes is given by 


cos^ 26 >k = 


2a‘^P 


(45) 

The analysis of these equations shows that nodal points 
appear in a set of 8 and that they are not located on the 
FSs of the normal state, although our numerical calcula¬ 
tions show that they remain very close to it. The loca¬ 
tion of the nodes with respect to the original FSs varies 
depending on the sign of a. If a > 0, the nodes are lo¬ 
cated inside of both unhybridized FSs but outside of the 
smaller FS. Instead, if a < 0 the nodes are outside the 
unhybridized FSs but inside the larger reconstructed FS. 

In both cases, the behavior is qualitatively the same 
as in the limiting case of circular pockets and is summa¬ 
rized in Fig. where we show the result of numerical 
calculations. Increasing A shifts the nodes toward the kx 


and ky axes. The critical value of A that causes the nodes 
to merge along these directions is enhanced by the ellip- 

ticity and is given by Ac = ^(A^ + {a‘^ — 1). At any 

larger A the nodes disappear. In this sense, the eccen¬ 
tricity of the pockets tries to prevent the disappearance 
of the nodes. 

Regarding the gap structure, we note that for large 
A the inter-pocket pairing term is irrelevant as the re¬ 
constructed FSs are far apart from each other. In this 
limit, the gap on the reconstructed pockets is given by 
A(1 ± /k^k), where |/k^k| ^ 1- Thus, the phase of the 
gap function is equal and uniform over the the recon¬ 
structed FSs. At smaller A, the phase of the gap along 
the FSs cannot be determined as the pairing involves 
both intra-pocket and inter-pocket terms. Judging from 
the large A limit, it seems natural to define the gap with 
equal sign on both FSs also at intermediate A, see Fig. 

m 


C. Inter-pocket pairing and hopping (/3 / 0, A / 0) 

This case is the most generic one. Like before, we 
switch to the reconstructed a and b pockets after hy- 
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Small inter-pocket 
hopping and interaction 




Figure 9: Gap structure for elliptical pockets with inter-pocket hopping and pairing. This case incorporates all the features 
seen in Figs. |^and[^ The main difference is that the condition for additional nodes does not depend solely on the ellipticity, 
but also on the hopping parameter and the sign of the inter-pocket pairing. 


bridization. The pairing terms in the Hamiltonian can 
be rewritten in terms of a and b operators and take the 
form 


^A+/3 — ^ ^ 


o ^ ^a6(k) <^kcK^Lk/3 + ^ka^- 


k/3 


Id 


a/3 


^^a/3 

H.C., 

( 46 ) 


The dispersions are given by where 




g£)" + kbG2A2 ^1 + ^2) +2/?2 


( 50 ) 


((eb'-(^k)V 4 AVk 2 /k) 

- 2 /?A 5 k (?£ - eb + (eg - eb /k^k] 


+/32 


4A2 + (eg-ebTi 


( 51 ) 


where 


Aaa(k) = A(1 - ^k/k) - /^^k, 
AbbO^) = A(1 + ^k/k) + 
AabO^) — A^k^^k H“ /^/k- 


( 47 ) 

( 48 ) 

( 49 ) 


Once more, we search for nodes in the dispersion . 
The radial position of the nodes is determined by the 
condition 


(^^k^k — 0'A‘^) ((5k + Q^^k) 

+ ^ {a^A^ - 6l)+a{l3^ - X^)6k = 0. ( 52 ) 
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The solutions to this equation (C, are needed to find 
the angular position of the nodes: 


cos^ 26>k 


_ F{U) 

(S^ + ’ 


(53) 


where 


F(C, S) = (S^ - a^A^) -h/]^ -A^ - A^) 

+ 4aAS (/3A - A^ . (54) 


J 


Solving for ^ in terms of S, we find that the solutions to Eq. (52) are 




1 1/2 > 

(a^A^ +S^f+ 4as(^ - a^A^) ^ + aS (A" - /?") ^ 


(55) 


The angular location of the nodes for these solutions is given by 
1 


cos^ 26>k = 


2a^S^ 


-a^A^ -h2aS/]A/AT 


(a^A^ +Pf+ 4aS^ (P - PA^) ^ + aS (A^ - /3^) ^ 


1/2 > 


(56) 


In this general case the interplay of the different pa¬ 
rameters is considerably more complicated that in the 
previous limiting cases, but it does not produce any new 
features. We find that in general the position of the nodal 
points depends on the signs of both P and a, not only 
their magnitudes. The dependence on the sign of (3 comes 
from bilinear terms of the form while the dependence 
on the sign of a is a consequence of the ellipticity of the 
pockets. 

The general behavior of the nodal points is summarized 
in Fig. We found by numerical analysis that in general 
increasing A tends to shift the nodes toward the and 
ky axis as usual, where they merge and disappear at a 
critical value of A. In this case, the sign of the order 
parameter is the same on both FSs. Increasing \P\ instead 
shifts the nodes toward the diagonals kx = Aky. At |/3| = 
A there are two possibilities, as seen in the limiting case 
of A = 0. The first is that the eight nodal points merge in 
pairs at the diagonal lines, disappearing for |/3| > A. This 
happens at small S, i.e., for small eccentricity. The second 
is that the original nodal points do not reach the diagonal 
lines at this value of /d, but instead four new nodes appear 
at those lines. Increasing \P\ further causes these four 
new nodal points to split into pairs, and moves the old 
and new nodes toward each other. At a threshold value of 
P they merge and then disappear. This second scenario, 
with 16 nodal points at intermediate is realized at 
larger J, i.e., at larger eccentricity. In either case, the 
merging of nodes caused by large P means that the gap 
function has opposite signs on the electron pockets. 

One important difference with the case of A = 0 is 
that the condition for developing additional nodal points 


is more complicated since the sign of P and the value of 
A also play a role. It is clear from the equations that the 
additional nodes are more likely to develop for < 0, 
for 5 > a A and for small A. The exact conditions when 
additional nodes appear are given by rather involved for¬ 
mulas and we refrain from presenting them. 

V. CONCLUSIONS 

In this paper we have investigated the effect of hy¬ 
bridization of the two electron pockets on the gap struc¬ 
ture in FeSCs. We considered the case when the dom¬ 
inant pairing interaction is between hole and electron 
pockets and it yields an s~^ gap with accidental nodes 
on the electron pockets. Our goal was to understand how 
accidental nodes move once we include the hybridization. 
We argued that for an s~^ superconductor hybridization 
gives rise to two effects - hopping between electron FSs 
and the appearance of an additional pairing term which 
describes inter-pocket pairing. Each of these two effects 
shifts the position of the nodes and at large enough hy¬ 
bridization the nodes eventually disappear. However, the 
evolution of the nodes and the gap structure of the re¬ 
sulting no-nodal state is different, depending on whether 
the inter-pocket hopping or the inter-pocket interaction 
is stronger. In the first case, the resulting state has the 
same sign of the gap on both reconstructed FSs. In the 
second case, there is a sign change of the superconducting 
gap between the inner and outer FSs. We also found that 
the evolution of the nodes with increasing inter-pocket 
pairing interaction is rather non-trivial, and in the inter- 
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mediate regime the number of nodal points may increase 
from 8 to 16. We also found that the eccentricity of the 
pockets enlarges these critical values, partially protecting 
the nodal points from disappearing. The bottom line of 
this analysis is that strong hybridization lifts accidental 
nodes, but the resulting superconducting state may be 
highly non-trivial, particularly when the dominant effect 


of hybridization is the emergence of inter-pocket pairing 
potential. 
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